Duality between quantum and classical dynamics for integrable billiards 
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We establish a duality between the quantum wave vector spectrum and the eigenmodes of the 
classical Liouvillian dynamics for integrable billiards. Signatures of the classical eigenmodes appear 
as peaks in the correlation function of the quantum wave vector spectrum. A semiclassical derivation 
and numerical calculations are presented in support of the results. These classical eigenmodes can 
be observed in physical experiments through the auto-correlation of the transmission coefficient 
of waves in quantum billiards. Exact classical trace formulas of the resolvent are derived for the 
rectangle, equilateral triangle, and circle billiards. We also establish a correspondence between the 
classical periodic orbit length spectrum and the quantum spectrum for integrable polygonal billiards. 
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Two sets of spectra can be associated with any dynam- 
ical system. The quantum spectrum, represented by the 
eigenvalues of the Schrodinger equation, is a set of real 
numbers for a closed system, while for an open system, 
the eigenvalues are complex and are called resonances. 
The classical dynamics can be described in terms of the 
spectrum of eigenvalues of the Liouville operator. For 
hyperbolic systems, the classical spectrum is comprised 
of the so-called Ruelle-Pollicott resonances [1, 2] which 
determine the time evolution and relaxation of classical 
correlations. 

Since the birth of quantum mechanics the correspon- 
dence between the classical and quantum properties has 
been examined from different perspectives. While the 
classical dynamics in the phase space is governed by the 
Liouville equation, the corresponding quantum dynam- 
ics in phase space is governed by the Moyal equation of 
the Wigner function. In the classical limit, the higher 
order terms of H in the Moyal equation vanish, and one 
retrieves the classical Liouvillian dynamics. This corre- 
spondence in the phase space properties has been studied 
by Brumer and collaborators [3]. In this paper, we obtain 
three principal results that establish a duality between 
quantum and classical dynamics in integrable billiards. 

For billiard systems, it is natural to work in wave vec- 
tor space [4, 5] . We find that the two-level correlation of 
billiard systems in wave vector k space is invariant with 
respect to the correlation interval. Peaks are observed 
in the quantum spectral correlations which are shown to 
be centered at the classical eigenmodes of the Liouvil- 
lian dynamics. Thus we establish a duality between the 
quantum wave vector spectrum {k n } and the classical 
eigenmodes {7™} for integrable closed billiard systems, 



C[f({k n })} = g({ 7n }). 



(1) 
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Here f{{k n }) and g{{"f n }) are certain spectral functions 
with C[f] an appropriately defined quantum correlation. 
The quantum wave vector spectral correlation directly 
leads to the classical eigenmodes. The results are numer- 
ically demonstrated for example systems including the 
two-dimensional rectangle, equilateral, and the circle bil- 
liards. A semiclassical derivation is provided supporting 
the results. Generalization to higher dimensional inte- 
grable systems is straightforward. 

It is well known that the quantum spectrum can be 
calculated from classical periodic orbits (POs) through 
the Gutzwiller trace formula. Here we obtain an exact 
inverse result, viz., that the classical PO length spectrum 
can also be expressed in terms of the quantum wave vec- 
tor spectrum, for integrable polygonal billiards. Finally, 
we also obtain exact classical trace formulas for integrable 
billiards. 

In Sec. II, the quantum correlation is defined and the 
billiard systems under consideration are introduced. The 
corresponding classical dynamics are discussed, the clas- 
sical eigenmodes and exact classical trace formula are 
obtained in Sec. III. A semiclassical derivation of the 
quantum correlation is provided in Sec. IV. Detailed 
numerical calculations supporting the main results are 
presented in Sec. V. Exact relation between the classi- 
cal POs spectrum and the quantum spectrum is estab- 
lished for integrable billiards in Sec. VI. Finally, Sec. VII 
presents a summary and some remarks on the results. 



II. QUANTUM EIGENVALUES AND 
CORRELATIONS 



In this paper, we consider the wave vector density of 
states (DOS) as our quantum spectral function which is a 
sum of ((-functions p(k) = J2 n S(k—k n ). In order to study 
the quantum spectral correlation, we use a Lorentzian- 
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smoothed DOS 
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Here e is a small width. In the limit e — > 0, one gets 
the stick spectrum. The continuous part of the DOS 
for two-dimensional closed billiard systems is obtained 
from the Wcyl law as (p(k)) = (A/27r)fc ± L/4ir with A 
and L the area and perimeter, respectively. The plus 
sign is for Neumann boundary condition (NBC) while 
minus sign for Dirichlet boundary condition (DBC). The 
fluctuation part of DOS is 6 p e (k) — p e (k)— (p(k)) . Define 
the following correlation 

C Pe (s) = {5pe(k)Sp € (k + s)) k . (3) 

Here 5p e (k) is with a window [K — A, K + A] with 
certain Kq and A such that A>(5, the mean level spac- 
ing. The range of average over k is from —A — min(0, s) 
to A — max(0, s). The range of s can be safely put 
as — A < s < A. The average over k is defined as 
(f{k)} = (b — a) -1 f(k)dk. Cross-correlation can be 
defined similarly for any two intervals \K\ — A, K\ + A] 
and [K2 — A, K2 + A] with arbitrary K\ and K 2 [6]. 

We remark that the correlation C Pe (s) in Eq. (3) is es- 
sentially the two- level correlation B,2(s). The difference 
is that here, there is no unfolding and no rcscaling by the 
mean level spacing of the energy levels. For billiard sys- 
tems, the correlation C Pc (s) in Eq. (3) is invariant with 
the choice of K n and A as long as one has enough quan- 
tum eigenvalues in the interval. For the Ricmann zeros 
[6-10], since the mean density is (p{k)} = (l/27r) ln(fc/27r) 
which is very flat and practically a constant for large k, 
there is no real difference between unfolding and not- 
unfolding. But for billiard systems, not-unfolding is es- 
sential to uncover the structures in the quantum spectral 
correlations. 

Here we focus on the quantum spectra of three intc- 
grable systems, the rectangle, equilateral triangle, and 
circle billiards. 

Rectangle billiard For rectangle billiard with sides a 
and b, the eigen wave vector is 

k mn = TrVWa) 2 + injby. (4) 

Here for DBC m, n > 1 and to, n > for NBC. 

Equilateral triangle billiard The equilateral triangle 
billiard is intcgrable only if all the sides have the same 
boundary condition. The eigenvalues are given by 

kmn = (47r/3a) \/to 2 + n 2 — mn (5) 

with a the side length. For the equilateral triangle with 
DBC, to > 2n > 2 while with NBC, to > 2n > 0. 

Circle billiard The circle billiard is integrable. The 
spectrum {k mn } is given by 

J m (k mn a) = 0, DBC, 

J' m (k mn a) = 0, NBC. (6) 

Here a is the radius of the circle and J m (x) is the Besscl 
function of the first kind. 



III. CLASSICAL EIGENMODES AND TRACE 
FORMULA IN LIOUVILLIAN DYNAMICS 

We first discuss the classical eigenmodes of the Liou- 
villian dynamics. For Hamiltonian systems, the phase 
space density distribution g(qi,Pi,t) is governed by the 
Liouville equation 

d 

■^Q{qi,Pi,t) = Lg(qi,pi,t) (7) 

where L = {H, } is the Liouville operator. Since we are 
interested in the spectrum of L instead of the detailed 
solutions of the equations of motion, it is better to switch 
to the action-angle variables which are related to qi and 
Pi through a canonical transformation [11]. In the action- 
angle variables, the Hamiltonian is a function of actions 
U = (27r)~ 1 §pidqi, H = H{Ii). For fixed Ii, the phase 
space is a N-dimensional torus with area (2^)^ for a N- 
dimensional integrable system. The equations of motion 
for the angle variables $i are dfii/dt = u>i with LOi = 
dH / dli, which give solutions as linear functions of time. 
These frequencies have dependence on action variables Ii 
and are in general continuous. They are often associated 
with certain trajectories. Since the Hamiltonian is cyclic 
in the Liouville operator L assumes a much simpler 
form in terms of action-angle variables 

- x - dH d \ - d 

The eigenvalues of operators ig§- all take integer val- 
ues for motions on the torus. Thus u>i are the primitive 
eigen frequencies of the Liouville operator L on the N- 
dimcnsional torus. 

Except for one-dimensional systems, these frequencies 
u>i are not the eigenmodes of the Liouville operator on 
the energy surface which has 2N — 1 dimensions. In or- 
der to obtain the classical eigenmodes of the system, an 
ensemble average must be performed over the whole en- 
ergy surface. A good way is to look at the trace of the 
classical evolution operator 

tr e Lt = - ^2 J dpexp(in ■ cot). (9) 

Here dp is the Liouville measure on the energy surface 
in phase space [12]. For two-dimensional integrable sys- 
tems, one has 

dp = 4ir 2 S(E - H{h,I 2 ))dhdI 2 . (10) 

Here the factor 47r 2 is from the integration over the two 
angle variables since any integrand we will consider for an 
integrable system is independent of the angle variables. 
The total measure is p = Air 2 J J dIidI 2 S(E — H) = 
dV/dE with V the phase space volume. The Laplace 
transform will give the trace of the resolvent 

/•OO 

g(z) = tr(z - Ly 1 = / tie Lt e- zl dl. (11) 
^0 
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Here for billiards, the particle speed v is set to be 
unity, thus time t is identified with length / through- 
out the paper. The classical eigenmodes will show up 
as singular points in the trace of the resolvent. For 
hyperbolic systems, the trace of the resolvent has sim- 
ple poles which are called the Ruelle-Pollicott resonances 
g(z) = Y,n d ™/( z ~ 7n) with d n the degeneracy. For 
generic systems, the classical trace g(z) will have various 
kinds of singularities which can be associated with the 
classical eigenmodes. 

Rectangle billiard For a particle with mass m mov- 
ing inside a rectangle billiard, the two actions are I x — 
pa cos a /it and I y = pb sina/n with p = \/2mE and 
a the angle of the velocity with respect to the x-axis. 
One thus gets u> x — (irv/a)cosa and u> y = (7rw/6)sina 
with v the particle speed in the billiard. Since the Jaco- 
bian determinant \d(I x ,I y )/d(p, a)\ is pab/ir 2 , so one has 
dfi = Aab5(E — p 2 /2m)pdpda. For the rectangle billiard, 
the two actions should be all positive, so a e [0, n/2] and 
the total measure is yU = 2nmab. One obtains 
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Here k mn is given by Eq. (4). This result was also ob- 
tained by Biswas [13]. The Laplace transform of the 
Bessel function Jq(x) gives 
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(13) 



The classical eigenmodes are given by y mn — ik mn . 
These classical eigenmodes are discrete. The above trace 
formula of g roc (z) diverges as 2abk;y/ir with Un the cutoff 
nf h 

Equilateral triangle billiard For a particle in an equilat- 
eral triangle billiard, the classical dynamics is integrable. 
The phase space surface is a regular hexagon with paral- 
lel sides identified, and thus has the topology of a torus 
[14]. The classical eigenmodes are given by j mn = ik mn 
with k mn given by Eq. (5) and m,n = 0, ±1, ±2, • • •. The 
trace of the resolvent is 



I 00 
9tn{z) = - 



1 



3 m ,^-oo \A 2 + k2 m 



(14) 



The trace of the resolvent can be derived from the fol- 
lowing expression 



tr e 



Lt 



^o />: 

m,n— — oo 
^ oo 



i(47r/3a)[m cos a+n cos(7r/3— ot)]l 



da 



(15) 



Though this is similar to that for the rectangle billiard, 
we are unable to give a derivation nor a physical inter- 
pretation for this expression. 

Circle billiard For a particle in a circle billiard of radius 
a, the Hamiltonian is H = p 2 ./2m+p 2 ,/2mr 2 . The actions 
for the angular and radial freedoms are Ig = pg and = 
tt^ 1 (p 2 a 2 - Ig) 1 / 2 - ir^ 1 Ie &rccos(Ig /pa) withp= \/2mE 
[15]. We parameterize the two actions as 1$ = pa cos a 
and I r — pa (sin a — a cos a) /tt with < a < tt since I r is 
positive and Ig can take positive or negative value [12]. 
The frequencies are uo r — -mi/asma and ivg = av/a sin a. 
They are both positive as required. Since the Jacobian 
determinant is \d(Ig,I r )/d(p,a)\ = \pa 2 sin 2 a, so d[i = 
4na 2 d(E - p 2 /2m)ps'm 2 adpda and \i = 2ir 2 ma 2 . Thus 
one has 
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Scir(z) 



e t( ma +nK)l/a sin a ^2 ^ ^gj 



m.n— — oo 



S(z; m, n) 



with 
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E(z;m,n) = - - . 2 
tt J a^^ 1 " 



az sin 4 ada 



z 2 sin a + (ma + nir) 2 



(17) 



(18) 



Here z = 7 — is with 7 > 0. 

Simple expressions can be obtained for E(z; to, n) with 
m = 0. One has S(z;0,0) = 1/az. For n/fl, one has 
S(z; 0, n) = E(z/n; 0, l)/n with 



H(*;0,1) = 



2n + \Ja 2 z 2 + tt 2 az 



[tt + \/a 2 z 2 + tt 2 ] 2 Va 2 z 2 + tt 2 ' 
Here we have used the following integral 



(19) 



cos 2nada 



/ z \ 2 « 



1 + z 2 sin 2 a \f\ + z 2 \ 1 + Vl + z 



From the above expressions of S(z;0,n), one gets the 
classical eigenmodes 



7o„ = irnr/a, n > 0. 



(20) 



For to ^ 0, no simple expression is found for S(z; to, n), 
though the classical eigenmodes can still be obtained 
from its integral representation. For m,n / 0, the de- 
nominator in the integral representation of S(z;m,n) 
may have two roots for certain z being purely imag- 
inary and a e (0, n). If these two roots are equal, 
S(z; to, n) will diverge. That is, z = ±if(ao) with 
f(a) = (ma + nix) /a sin a and its first derivative f'(a) = 
{m—(ma+nir) cot a] /a sin a vanishes at ao- The classical 
eigenmodes are thus given by 



m.n— — 00 



Iran = imy/l + u 2 /a 



(21) 
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FIG. 1: Function 3(7 - is; 1,0) with 7 = 0.001 and a = 1. 
The solid line is real part 5RH(7 — is; 1, 0) and the dashed line 
is the imaginary part SsE(7 — is; 1,0). 



with u the solution of the transcendental equation 
u = mr/m + arctanw. 



(22) 



For to ^ and n — 0, one has ^(z; to, 0) = 
S(z/m; l,0)/m. 2(7 — is; 1,0) is free of divergence even 
when 7 -> and SRH(-is; 1,0) = for s < 1. We use 3? 
for the real part and for the imaginary part through- 
out the paper. For s > 1, is; 1,0) > and has a 
maxmum at s ~ 1.1525. The function 2(7 — is; 1, 0) with 
a small 7 is evaluated numerically and plotted in Fig. 1. 
The resonance-looking hump in is; 1, 0) gives rise to 
continuous classical eigenmodes and may be the reason 
for the apparent nonzero Lyapunov exponent [16] in the 
circle billiard and exponential short-time decay in circle 
billiard with small holes [17]. We emphasize that these 
classical eigenmodes bear no resemblance to the quantum 
eigenmodes. The first few discrete classical eigenmodes 
are listed in Table I. 



TABLE I: A few low- lying quantum and classical eigenmodes 
of the circle billiard. 



k mn for NBC 




kmn for DBC 








2.404825558 


1.841183781 


3.141592654 


3.831705970 


3.054236928 


4.603338849 


5.135622302 


3.831705970 


5.943387741 


5.520078110 


4.201188941 


6.283185307 


6.380161896 


5.317553126 


7.221364982 


7.015586670 


5.331442774 


7.789705768 


7.588342435 


6.415616376 


8.460248568 


8.417244140 


6.706133194 


9.206677698 


8.653727913 


7.015586670 


9.424777961 


8.771483816 


7.501266145 


9.671732536 


9.761023130 


8.015236598 


10.566779006 


9.936109524 



For generic systems, the frequencies uji can not be un- 
entangled with each other since the phase space is not 



a torus. Eq. (9) is not useful to obtain the classical 
eigenmodes. However the trace of the resolvent can be 
expressed in terms of POs which are the periodic solu- 
tions of the classical equations of motion. For hyperbolic 
systems, Cvitanovic and Eckhardt [18] obtained the clas- 
sical trace formula for the resolvent which is similar to the 
Gutzwiller semiclassical trace formula [19]. The classical 
trace formula for polygonal billiards is derived recently 
by Biswas [13]. Nevertheless, exact classical trace formu- 
las can be derived directly from Eq. (9) for integrable 
billiards. 

For the rectangle billiard, the fluctuation part of the 
classical trace formula can be directly derived from Eq. 
(12) as shown in the Appendix as 



r / \ 2ab 

ogr CC {z) = — 

7T 



m,n— 1 



4 E 
+2E 



m— 1 



zL„ 



-zLn 



J m0 



J 0rn 



Here the length of the POs is given by 



ma) 



{nbf 



(23) 



(24) 



For the equilateral triangle billiard, along the same 
lines of derivation for the rectangle billiard, one gets 



<5fftri(z) 



V3a 2 
2tt 



-zLr, 



6E 



,-zimO 

LmO 



with 



(25) 



(26) 



For the circle billiard, the classical trace formula can 
be derived directly from Eq. (16) as shown in the Ap- 
pendix as 



Sgcir(z) 
with 



1 ^ y ^ y J inn 

n—1 m—2n 



sin 4 (mr /to) Ln 



(27) 



f mn = 2 for to > 2n, 
= 1 for to = 2n, 
L mn — 2ma sin(n7r/ro). 



(28) 



In summary the above represents exact trace formulas 
for the classical dynamics of integrable billiards. 



IV. SEMICLASSICAL DERIVATION OF 
QUANTUM CORRELATIONS 

A semiclassical consideration can lead to the under- 
standing of the quantum correlation. For large fc, the 
semiclassical DOS for integrable polygonal billiards is [20] 



S P e(k) 



9^E a P 



^{ik — e)Lp—i7r/A 



(29) 
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The auto-correlation C Pe (s) can be written as the sum 
of diagonal and off-diagonal terms, C Pc (s) = C^ iag (s) + 
C*g(s). For integrable systems, the actions are uncorre- 
lated since there is a direct relation between actions and 
the quantum spectrum. This leads to the cancelation of 
the off-diagonal terms, 



For the circle billiard, one has 



s)L p 



(30) 



Note that the above summation is over POs with dif- 
ferent L p . Different POs with the same L p are grouped 
together. For the rectangle billiard with irrational aspect 
ratio, ftp = Aab for all PO families L mn with nonzero m 
and n except for L m0 and L 0n . Thus 



C P Js) 



(fflp) K 
8tt 2 



%i6g rcc (2e - is). 



(31) 



With normalization at s = 0, C Pe (s) and $l5g(2e — is) 
are almost identical with each other though the relative 
heights of their peaks may be slightly different for rect- 
angle billiards with rational aspect ratio or the equilat- 
eral triangle billiard. This is due to the fact that a p is 
no longer the same for all POs with different length for 
these billiards. This gives rise to extra fine structure in 
C Pt (s). We point out that the appearance of fine struc- 
tures in C Pe (s) is not an indication of discrete nature of 
the quantum spectrum as claimed in [21, 22], rather that 
a p is not uniform. 

For the circle billiard, there is no exact trace formula 
for p(k) due to the lack of analytic expression for the 
zeros of the Bessel functions. For the billiard with DBC, 
the fluctuation part of the semiclassical DOS is [20, 23] 



n— 1 m—2n 



(32) 

with 3> mn = (k + ie)L mn + v m , v m = (6m - and 
fmn given in Eq. (28). This trace formula will give the 
WKB quantization [15, 20] 



fcon = (« + 3/4)7T 



which approximates the zeros of Jq{z) and 



m 



VI + u 2 



(33) 



(34) 



with 



u = (n + 3/4)7t/to + arctanw 



for the approximate zeros of J m (z) with m > 0. For the 
billiard with NBC, the factor 3/4 should be replaced by 
1/4. With either boundary condition, the diagonal part 
of the auto-correlation is 



oo oo 



sin (rwr/m) (2e _ <a)zLn 

n—1 m—2n 



Cdiag( S ) 



47ra 2 K 
8tt 2 



$t5g c i T (2e — is). 



(36) 



(35) 



Thus the quantum correlation is determined by the clas- 
sical eigenmodes. 



V. NUMERICAL CALCULATION OF 
QUANTUM CORRELATIONS 

Numerical calculations of C Pe (s) were carried out to 
confirm the results discussed here. Results for the rect- 
angle billiard with DBC are shown in Fig. 2. For the 
equilateral triangle billiard of side length a = 1 with 
DBC, the spectral correlation is shown in Fig. 3. Since 
the spectrum of the rectangle and equilateral triangle bil- 
liards with NBC is almost the same as the counterparts 
with DBC, their spectral correlations are essentially the 
same. The spectral correlations for the circle billiard of 
unit radius is shown in Fig. 4. About 10 6 quantum 
eigenvalues are used in the above calculations. 

As seen from the above figures, the auto-correlations 
C Pt (s) are almost the same for different spectra segments, 
confirming the correlation invariance. It is quite obvious 
that the peaks of the quantum wave vector correlation are 
located at s = Qjmn with 7 m „ the classical eigenmodes. 

In Fig. 2, the rectangle billiard has an irrational aspect 
ratio. In this case, almost all the PO families have the 
same summation weight. With fixed quantum spectral 
width e, for a finite segment of quantum spectrum, the 
difference between the correlation C Pt (s) and the classi- 
cal trace $lSg(2e — is) can be viewed as just noise if both 
are normalized at s = 0. This noise will be reduced and 
eventually disappears if either the number of the quan- 
tum eigenvalue in the quantum spectral segment or the 
width e increase. However for the equilateral triangle bil- 
liard as shown in Fig. 3, the difference between C Pt (s) 
and 5Wg(2e — is) will persist even if the center K of the 
quantum spectral segment is pushed to infinity. Actu- 
ally, there are five curves with different Ko plotted in 
Fig. 3. They are hardly distinguishable from each other. 
The classical trace %i5g(2e — is) has no structure between 
classical eigenmodes. The presence of small peaks be- 
tween classical eigenmodes in C Pe (s) is what we called 
fine structure in the previous section and is due to the 
nonuniformity of a p . 

An interesting point about the circle billiard is that 
there are some resonance-like humps in the correlation as 
shown in Fig. 4 between the origin and the first nonzero 
discrete classical eigenmodes 701 = iir. The hump be- 
tween [l,7r/2] is very prominent. These humps may be 
the precursor of Ruelle-Pollicott resonances in chaotic 
systems such as the Bunimovich stadium. 

Two comments are in order. The first is that though it 
is tempting to make a one-one correspondence between 
the quantum and the classical eigenmodes for integrable 
systems, actually this is not true. For any segment of 
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FIG. 2: Auto-correlation C Pe (s) of the rectangular billiard with sides a — tt and b — \fn. Here e — 0.005. 10 6 quantum 
eigenmodes with < k mn < 1500 are included in 5p e (k). The auto-correlations C Pt {s) for Sp e (k) in four equal intervals [0, 375], 
[375, 750], [750, 1125], [1125, 1500] are calculated. They are normalized such that C Pe (0) = 1 and all collapse to the same curve. 
Vertical dashed lines indicate the location of the classical eigenmodes. 



the quantum spectra, the auto-correlation will lead to 
the same low-lying classical eigenmodes as exemplified 
by the invariance of the quantum spectral correlations. 
The second is that in general, there is no correspondence 
between the trace of the classical project operator and 
the quantum Neumann spectrum. One can clearly see 
this for the circle billiard. For the rectangle and equi- 
lateral triangle billiards, the fact that the quantum NBC 
spectrum is also the classical eigenmodes is a coincident 
instead of a rule for polygonal billiard. One can view 
the circle billiard as the limit of equilateral polygon with 
increasing number of sides. 



and classical dynamics such that the quantum DOS is 
expressed in terms of classical POs. As a complement 
to the Gutzwiller trace formula, the classical PO length 
spectrum can also be expressed in terms of quantum wave 
vector spectrum. This classical POs spectrum is crucial 
for the calculation of classical eigen spectrum using the 
Cvitanovic-Eckhardt classical trace formula [18]. 

For integrablc billiards, the classical trace formulas we 
derived in Sec. Ill assume the form 



p. e - ZL p 



(37) 



VI. CORRESPONDENCE BETWEEN 
CLASSICAL PERIODIC ORBIT LENGTH 
SPECTRUM AND QUANTUM WAVE VECTOR 
SPECTRUM 

In the previous sections, we established a duality be- 
tween quantum and classical eigen spectra through quan- 
tum correlations. This duality was also established for 
hyperbolic systems [24]. We believe this can also be ap- 
plied to generic systems. 

On the other hand, the celebrated Gutzwiller trace 
formula indicates a direct duality between the quantum 



with a p the area visited by the PO family with length 
L p . This formula is also valid for generic polygonal bil- 
liard [13]. In order to obtain the classical eigenmodes, the 
number of PO families defined as N(l) = J2 P a p Q(l — L p ) 
should be known. Here Q(x) is the step function. N(l) 
has the unit of billiard area. It was proved that the 
POs are dense for rational polygonal billiards [25, 26], 
N(l) oc I 2 . Recently Biswas obtained an expression of 
N(l) in terms of the quantum spectrum with NBC [27]. 
Actually, N(l) can also be expressed in terms of the quan- 
tum spectrum with DBC on the same basis. Here we give 
a different derivation of this connection and reveal the 
domain of its validity. Often the importance of isolated 
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Q- 0.1 




FIG. 3: Auto-corrclation C Pe (s) of the equilateral triangle billiard with side length a — 1. About 10 6 quantum eigenmodes 
with < k mn < 5390 are included in 8p € (k) with e = 0.059. The auto-correlations C Pe (s) for 5p e (k) in five equally divided 
intervals are calculated, normalized at the origin, and plotted together. Vertical dashed lines indicate the location of the 
classical eigenmodes. 



POs is neglected in this connection. We find that they 
arc key to determine the proliferation of the PO families. 

Since wc focus on integrable billiards, diffractive POs 
and non-periodic orbits [28-30] are absent, the wave vec- 
tor DOS of integrable polygonal billiards is 



p(k) = (p(k)) + c S(k) + — ^2a p J (kL p ) 



— bp cos kip. 

P 



(38) 



Here cq is the corner term and b p is the weight for isolated 
PO with length £ p [20]. For two-dimensional billiards 
with DBC, the average DOS is (p(k)} = Ak/2n - L/4tt. 
Since p(k) = ^2 n S(k — k n ), multiplying both sides of Eq. 
(38) by 2irlJ\(kl)/k and integrating over fc from zero to 
infinity, one gets 

T — R °° i 

Nil) = j(a-2 Co )l 2 + —^l-A + 2nlY / 7rMkJ) 



n— 1 p' 



(39) 



The summation over p' is for the primitive isolated PO 



and 



p' 



(40) 



p' 



In the above derivation, use has been made of the follow- 
ing identities [31] 



f dk.h{kLp).h{kl) = y9(/ - Lp) 
Jo ' 

L 



dk 1 

— cos(H p ) Mkl) = -j^Jl 2 - I 



3? ^2 \fx 2 - n 2 = —x 2 - 2 X + 2 E - J i( 2n7TX )- 

n=l 



n=l 



The terms wal 2 /4 and Pl/2 in N(l) are from the first two 
terms of the third identity. Thus one has the average 

(N(l)) = l( a -2c )l 2 + ^ll-A. (41) 

Hence the proliferation of PO families is determined by 
the isolated PO and the corner term in the quantum trace 
formula (38). This reveals an intimate connection be- 
tween them. 
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FIG. 4: Auto-correlation C Pc (s) of the circle billiard with radius a = 1. About 2.3 x 10 6 quantum eigenmodes with < 
kmn < 3000 are included in 8p e (k) with e = 0.005. The auto-correlations C Pe (s) for 5p e (k) in four equally divided intervals 
are calculated and normalized at the origin. Vertical dashed lines indicate the location of the classical eigenmodes. The two 
humps in the interval [1, 3] correspond to 5tH(2e — is; 1, 0) and 5?E(2e — is; 2, 0), respectively. 



Exact expression of N(l) can be obtained for the rect- 
angle and equilateral triangle billiards. For the rectangle 
billiard with DBC, the exact trace formula [20] will give 
c = 1/2, a = 2, ft = L = 2(a + b). One has exactly 

N{1) = jl 2 -ab + 2irl ^ — Ji(fc ro „Z) 

4 m,n=l mn 

+(V - laJUnnl/a) + bJUmrl/b)]. (42) 

n— 1 

Thus (-/V(O) = 7rZ 2 /4 - ab. We stress that there is no 
linear term and the constant term is negative in (N(l)), 
contrary to the expression obtained by Jain [32]. As a 
comparison, the number of quantum states is 



M(k) 



ab . 
4ir 



p ' 

1 00 1 

> — (sin 2nka + sin 2nkb) 

2tt ^ n 

n=l 



(43) 



which is obtained from Eq. (38) through an integration. 
Here the summation over p is given by Eq. (23). One 
has (Af(k)) = abk 2 /4w - (a + b)k/2ir + 1/4 as expected. 



For the equilateral triangle billiard with DBC [20], one 
has c = 2/3, a = 2, /3 = L = 3a. Thus 

m = ^_^ +2)ri ^j l(W) 



n=l 



3 °° a 
+-^lY] -Ji(4ri7r//3a) 

n— 1 



(44) 



with (N(l)) = 7r/ 2 /6 — V3a 2 /4. The number of quantum 
states is 



Af(k) 



1 . , 1 

^ — sin(3nfca/2). 



(45) 



Here the summation over p is given by Eq. (25). One 
has (JV(fc)} = V3a 2 /c 2 /167T - 3afc/47r + 1/3. 
Similarly, one gets 



p(i) ee x: - i P ) 
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= |0-2co)- Ak N ,h{k N l)/l + ^^- 

+2n MkJ) + ^ E E f;M2mrl/£ p/ ), 



n=l 



n=l p' 



(46) 



S(l) = J2 a p S ( l - L p) 



|(a - 2c )l - AkNJ^kNl) + 



N 



+2irl J2 -MM + ^ E E j^Jo^l/ip'). 

n—1 n—1 p' P 

(47) 

Here fc^v is a large cutoff in the quantum spectrum. Ex- 
act expressions of P(l) and S(l) can be obtained for the 
rectangle and equilateral triangle billiards. 

For billiard with uniform boundary condition, a p in the 
quantum trace formula (38) are all positive, the quantity 
P(l) is just the classical quantity (7r/2)tr exp(— Lt). The 
Laplace transform of P{1) will give the fluctuation part 
of the classical trace g(z) 



Sg(z) = (a- 2c )/z - 



N 

E 



2A 



z 2 + k% 



EE 



2y/y 



(48) 

Here we ignored the term (L — (3)/2l in P(l) since it 
vanishes for integrable polygonal billiards and made use 
of the Laplace transforms of Bessel functions 

J r x Jr(k N l)e- zl dl = ^-(yW^-z), 
f°° 1 

L Mk -' ,er ""' - T^m 

Thus the classical eigenmodes are 7„ = ik n and 7„ = 
i2mr/£ p > with £ p i the primitive isolated PO without rep- 
etition. 

The results for N (I) , P(l), S(l) establish a correspon- 
dence between the classical PO length spectrum and the 
quantum spectrum for integrable polygonal billiards. 

For the circle billiard, no exact expression of N(l) can 
be obtained since we don't have an exact quantum trace 
formula. Nevertheless, an approximation of N(l) could 
be derived similarly through the semiclassical trace for- 
mula. 

We remark that the above equations may be applicable 
to non-integrablc billiards only if the contribution to the 
quantum trace formula from diffractive POs and non- 
periodic orbits can be safely ignored. This is not the case 
for generic polygonal billiards. For example for the 7r/3 
rhombus billiard, diffractive POs are nearly as dense as 



the PO families. The classical spectrum is thus different 
from the quantum spectrum though they may partially 
overlap with each other. 



VII. CONCLUSIONS AND REMARKS 

In this paper, we directly obtained the classical eigen- 
modes of the Liouvillian dynamics and exact classical 
trace formula of the resolvent for integrable closed sys- 
tems including the rectangle, equilateral triangle and the 
circle billiards, and showed that the peaks in the quan- 
tum wave vector spectral correlation coincide with the 
classical eigenmodes. We also established a correspon- 
dence between the classical PO length spectrum and the 
quantum wave vector spectrum for polygonal integrable 
billiards. 

We have shown earlier in open chaotic systems, specif- 
ically the n-disk billiards, that the auto- correlation of 
the resonant quantum spectrum carries the fingerprints 
of the classical resonances, also known as the Ruelle- 
Pollicott resonances. For the Ricmann zeros which are 
the quantum eigenmodes on negative constant curvature, 
the spectral correlation also leads to the "classical" eigen- 
modes [6, 10]. We speculate that this also holds for 
chaotic and generic billiard systems. A direct duality 
between classical POs spectra and quantum spectra is 
also obtained for integrable billiards. 

While for a quantum system, there are many ways to 
solve the Schrodingcr equation, the classical spectrum is 
very difficult to get at. Direct numerical simulation of 
tracing classical trajectories is of limited use due to the 
presence of the Lyapunov exponent and computer round 
off error. The trace formula in terms of classical POs pro- 
vide another correspondence between the quantum and 
classical dynamics. The usefulness of these trace for- 
mula to obtain classical [18] and quantum spectra [19] 
is marred by the divergence due to the exponential pro- 
liferation of POs. Though the classical eigenmodes are 
difficult to calculate for generic systems, our approach 
provides a simple way to obtain them. 

Our results can be readily tested in wave-mechanical 
experiments, such as the microwave analogs of quan- 
tum billiards [33]. There the measured transmission is 
T(k) = J2 n c n/[{k - k n ) 2 + el]. Earlier work has shown 
in n-disk open systems that the peaks of the correlation 
Ct{s) = (T(k)T(k + s)) k are indeed located at s — ^s-f n 
[24]. In another words, the classical eigenmodes were ob- 
tained from the experimentally measured quantum spec- 
trum for this open system. Similar experiments on closed 
billiards should yield results that can be compared with 
those discussed in the present paper. We have examined 
simulations of the experimental data by taking c„ as ran- 
dom numbers in [0, 2e/ir] and set e„ = e, the same width 
for all eigenvalues as we used in Sp e (k). We find that 
spectral correlations are almost identical if c„ is constant 
(as in Sp e (k)) or uniformly distributed in [0, 2e/ir]. 

Classical eigenmodes play an important role in quan- 
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turn dynamics. For example they govern the time evolu- 
tion of a wave packet and are related to the problem of 
decoherence [34]. Consider the quantity 



A{t) 



ihk 2 n t/2m 



with m the particle mass, and a n is significant only for 
K - A < k n < K +A with K > A. This quantity is re- 
lated to the wave packet revival [35] . For two-dimensional 
billiards, the Heinsenberg time is tn = mA/h. For the 
wave packet moving with a speed v = hK /m, we de- 
fine the Heinsenberg length = vtu = (A/2n)K n . So 
one has Ih — ► oo if A — > oo (open systems) or Kq — > oo 
(short wavelength). One has Hk 2 t/2m = k n l — K l/2 + 



{A/4ir)(k n - K Q ) 2 (l/l H ). Here I = vt. 

last term can be ignored, thus |^4(£)| 2 
In the case a n = 1, one has 



For . 
E, 



< Z 



H 



a„.e 



the 

2 

ik n l 



\A{t)\- 



s)e isi ds ~ & e Lt . 



So before the Hcisenberg time, |^4(i)| 2 is governed di- 
rectly by the classical eigenmodes {7„} of the Liouville 
operator L. 

Boundary conditions play an important role on the 
dynamics. This was addressed recently by Biswas [36]. 
Non-identical boundary condition on adjacent edges of 
the polygonal billiards can lead not only to quantum 
splitting but also to quantum annihilation of PO families 
in the quantum trace formula. Such as in the case of a 
rectangle billiard with rational aspect ratio and only one 
side is with NBC. The most dramatic case is a square bil- 
liard with DBC on all sides except one. In this case, the 
quantum annihilation leads to a massive disappearance 
of PO families in the quantum trace formula. Similar to 
the case of quantum splitting, the spectral rigidity will 
also deviate from the Poissonian value due to quantum 
annihilation. 

Previous studies of the quantum properties of billiard 
systems were focused on the universality of level statis- 
tics, such as nearest level spacing, spectral rigidity, etc., 
and the quantum-classical correspondence was typically 
addressed in terms of POs [19, 37]. Here we take a differ- 
ent approach by examining correspondence between the 
classical and quantum spectra. In general, the classical 
eigenmodes arc different from the quantum eigenmodes. 
Except for some special systems [27, 38, 39], there is 
no direct correspondence between the two spectra. For 
systems in which all the POs have the same Lyapunov 
exponent, there is a self-duality between quantum eigen- 
values and the classical eigenmodes, such as the Riemann 
zeros. This self-duality leads to the resurgence [6-10] and 
prophecy [6] of quantum eigenmodes in the spectral corre- 
lation. For systems such as the rectangle and equilateral 
triangle billiards we considered in this paper, there is an 
approximate self-duality between these two spectra, the 
quantum eigenmodes are also the classical eigenmodes. 
Prophecy [6] are also observed in these systems. 



A general theory connecting the quantum fluctuation 
and the classical spectrum of the Perron- Frobenius opera- 
tor has been developed [40-42] for diffusive systems. Due 
to the difficulty of getting the classical eigenmodes for 
generic systems, the focus was shifted onto the connec- 
tion with random matrix theory. Since the spectral form 
factor is not self-averaging for clear systems, smoothing 
or ensemble average is required for the form factor of a 
single system to approach the prediction of random ma- 
trix theory [43]. Our previous study confirmed that for 
hyperbolic open systems, the quantum correlations arc 
determined by the classical eigenmodes [24]. In this pa- 
per, we have confirmed this for integrable closed billiards. 
Thus one is able to extend the so-called AAA-BK theory 
[40, 41] to a finite spectrum of a single closed system. 
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APPENDIX A: DERIVATIONS OF CLASSICAL 
TRACE FORMULAS 

In this appendix, we derive the trace formulas Eq. (23) 
and (27) from Eq. (12) and (16), respectively. 
For the rectangle billiard, one has explicitly 



tr e 



r, nTl/2 oo 
Lt / \ ^ ^im(iT /a)l cos a+in(7r /b)l sin a 

~ T J ^ 

u m,n— — oo 

Using the identity 

oo oo 

e 27Ttmx = S(x-m), (Al) 

m— — OO 

one gets for I > 



m— — oo 



8ab 



7T ' * 



(A2) 



m,n— — oo 



with 



/■tt/2 

F; ec „(I)= / 5(lcosa-2ma)6(lsma-2nb)da. (A3) 
Jo 

The above defined function F™ c n (l) is nonzero only if 
both m and n are non-negative. So we only need to 
consider to, n > 0. Using the property of the (5-function 
in the polar coordinate system [44] 

<5 2 (r-r ) = S(x - x )S(y - y ) 

= -S(r - r o )S(0 - 6 ) (A4) 
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we get for to, n > 0, 



1 



1 



r/2 



<5(^ — L mn ) 



+5(a - 7r + a m ) 

<5(a — a mn )da w i tri = arcsin(^/2m.a) for to > Z/2a, we get 

; 3 



(A9) 



(A5) F™ n (l) = 



with L mn — 2\Jm 2 a 2 + n 2 b 2 and a mn — tan -1 (rib / 'ma) . 
For F^ (l) with to > 0, since a m o = and J °° 5(x)dx = 
1/2, one gets F™ c (l) = 5(1 - L m0 )/2L m0 . Similarly one 
has F^(l) = 5(1- L 0n )/2L 0n and F^(l) = 5(1) /2l. So 
one gets 



<5(Z - Imp) £(l - L 0m ) 



m—l 



J m0 



J 0m 



+4E 



<5(^ — L mn ) 



(A6) 



m,n— 1 



The Laplace transform of the above expression without 
the first term 5(1) /I will give Eq. (23). 

For the circle billiard, the trace of the Liouvillian dy- 
namics is 



tre Zt = 8a 2 £ F% n (l) 



(A7) 



with 

pcir 



n (l) = f da sin 2 a5(—^- 2ma)5(-^- 2mra). 

Jo 



sin a 



sin a 



(A8) 

The above defined F™ T n (l) will be nonzero only for m,n> 
0. For to > 0, using 



5(- 



1 



sin a 



— 2ma) = 



2ma\ / '4m 2 a 2 — I 2 



5 (a - a m ) 



(2ma) 3 ^Am 2 a 2 - I 2 
I 3 

(2too) V4m 2 a 2 - I 2 



5(a m lsina m — 2mra) 
5(a m - nir/m) (A10) 



which is valid only for to > 2n. Since a m — nn/m will 
give I = L mn = 2masm(mr/m), thus 



5(a m - rm/m) = y/ Am 2 a 2 - L 2 mn 5(l - L mn ), (All) 
we obtain 



F m ,n(l) 



sin 4 (n7r/m) 



5(1 — L mn ). 



(A12) 



One can verify that this expression is valid for all to > 
n > 0. For to = and n > 0, one has fon(0 = 
while F^ (Z) = l 3 5(l)/ir(2ma) 4 for to > 0. If only the 
contributions from to, n > arc included, one gets 



tr e 



LI 



sin (mr/m) 



5(1 — Z/ mn ) 



n— 1 m— n 

CO CO 



n— 1 m-2n 



with 



a m n = 47ra 2 / m „ sin 4 (n7r/TO) 



(A13) 



(A14) 



and fmn given in Eq. (28). A Laplace transform will give 
Eq. (27). 
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